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LIOUVILLE-TYPE THEOREM FOR A FOURTH 
ORDER SUPERCRITICAL PROBLEM 
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Abstract. We consider Liouville-type and partial regularity re- 
sults for the nonlinear fourth-order problem 

A 2 u = \u\ p ~\ in R n , 

where p > 1 and n > 1. We give a complete classification of 
stable and finite Morse index solutions (whether positive or sign 
changing), in the full exponent range. We also compute an upper 
bound of the Hausdorff dimension of the singular set of extremal 
solutions. Our approach is motivated by Fleming's tangent cone 
analysis technique for minimal surfaces and Federer's dimension 
reduction principle in partial regularity theory. A key tool is the 
monotonicity formula for biharmonic equations. 



1. Introduction 

We study the following model biharmonic superlinear elliptic equa- 
tion 

A 2 u = \u\ p ~ 1 u inQ, (1.1) 

where Q C M. n is a smoothly bounded domain or the entire space and 
p > 1 is a real number. Inspired by the tangent cone analysis in 
minimal surface theory, more precisely Fleming's key observation that 
the existence of an entire nonplaner minimal graph implies that of a 
singular area-minimizing cone (see his work on the Bernstein theorem 
[12]), we derive a monotonicity formula for solutions of (II. ip to reduce 
the non-existence of nontrivial entire solutions for the problem (II. ip . 
to that of nontrivial homogeneous solutions. Through this approach 
we give a complete classification of stable solutions and those of finite 
Morse index, whether positive or sign changing, when Q = ~R n is the 
whole euclidean space. This in turn enables us to obtain partial regu- 
larity as well as an estimate of the Hausdorff dimension of the singular 
set of the extremal solutions in bounded domains. 



Key words and phrases. Monotonicity formula, stable or finite Morse index equa- 
tions, biharmonic equations, partial regularity. 
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Let us first describe the monotonicity formula. Equation fll.ip has 
two important features. It is variational, with energy functional given 
by 



it 



ip+i 



V 



and it is invariant under the scaling transformation 

, 4 

u (x) = \p~ 1 u(\x). 
This suggests that the variations of the rescaled energy 

-« f ±( A u) 2 - —\u\ p+1 



42+1. 
7" p-1 



with respect to the scaling parameter r are meaningful. Augmented 
by the appropriate boundary terms, the above quantity is in fact non- 
mcreasme;. More precisely, take u e W^iVL) n Lf^(n), fixiGfl, let 
< r < R be such that BJx) C B R (x) C f2, and define 
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where derivatives are taken in the sense of distributions. Then, we have 
the following monotonicity formula. 



Theorem 1.1. Assume that 

n > 5, 



p> 



n + 4 



1.2] 



n 



Let u G W lm (f2)flLf oc (fi) 6e a u>ea£; solution of (II. ip . TTien, £'(r; x, tt) 
non-decreasing inr e (0, i?) . Furthermore there is a constant c(n, p) > 
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such that 

^E(r;0,u)>c(n,p)r- n+2+ ^ J ^_i_ r -i tt+ ^ . (1. 3 ) 

Remark 1.2. Monotonicity formulae have a long history that we will 
not describe here. Let us simply mention two earlier results that seem 
closest to our findings: the formula of Pacard [23] for the classical 
Lane-Emden equation and the one of Chang, Wang and Yang [2] for 
biharmonic maps. 

Consider again equation (11. ip in the case where Q = R n , i.e., 

A 2 u = \u\ p ~ l u in R n . (1.4) 



Let 



+oo if n < 4 
Ps{n) = { n + 4 

if n > 5 



n - 4 

denote the Sobolev exponent. When 1 < p < ps(n), all positive solu- 
tions to (jl.4p are classified: if p < Ps{n), then u = 0; if p = ps(n), then 
all solutions can be written in the form u = Cn( \2 + \x- Xo \ 2 )~ T ~ ^ or some 
c n > 0, A > 0, x G M. n , see the work of Xu and one of the authors [33] . 
However, there can be many sign-changing solutions to the equation 
(see the work by Guo, Li and one of the authors [TB] for the critical 
case p — ps(n)). 

Here, we allow u to be sign-changing and p to be supercritical. In- 
stead, we restrict the analysis to solutions having finite Morse index. 
More precisely, define the quadratic form 

A u (0) := f \A(f)\ 2 dx-p [ lu^tfdx, V G H 2 (R n ). (1.5) 

JM. n JM. n 

A solution to (11. 4p is said to be stable if 

\A(j)\ 2 dx -p I H P ~V 2 <fc > 0, V <p G H 2 (W l ). 



More generally, the Morse index of a solution is defined as the maximal 
dimension of all subspaces E of H 2 (M. n ) such that A u ((f>) < in E\ {0}. 
Clearly, a solution is stable if and only if its Morse index is equal to 
zero. It is also standard knowledge that if a solution to (11.41) has finite 
Morse index, then there is a compact set /C C W 1 such that 

f \A(f)\ 2 dx -pi \u\ v ~ x <\) 2 dx > 0, V G if 2 (M n \/C). 
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Recall that if 

7 = : r, = 7(7 + 2) (7 - n + 4) (7 - n + 2), (1.6) 

p — 1 

then 

u s ( r ) = ^V(p-i) r -4/(p-i) ( L7 ) 

is a singular solution to (11.4p in IR n \ {0}. By the Hardy-Rellich in- 
equality with best constant [28] 

\A<l>\ 2 dx> n2{n ~ A)2 [ ^dx, ^eH\R n ), 
the singular solution u s is stable if and only if 

lb 

Solving the corresponding quartic equation, f ll.Sp holds if and only if 
P > Pc(^) where p c (n) > ps(n) is the fourth-order Joseph-Lundgren 
exponent computed by Gazzola and Grunau [14] : 

+00 if n < 12 



Pc(rc) = < ra + 2 - V« 2 + 4 - nVn 2 - 8n + 32 „ 

^^^=^^^=^^^^^^^= if n > 13 

n — 6— v^ 2 + 4 — n\/n 2 — 8n + 32 

Equivalently, for fixed p > ps{n), define n p to be the smallest dimension 
such that fll.8p holds. Then, 

(11. 8p -<=>- p > pdji) -<=>■ n > rip. 

The existence, uniqueness and stability of regular radial positive so- 
lutions to (11.41) is by now well understood (see the works of Gazzola- 
Grunau, of Guo and one of the authors, and of Karageorgis [141 [T71 12"0] ) : 
for each a > there exists a unique entire radial positive solution 
« a (|x|) to (11.41) with u a (0) = a. This radial positive solution is stable 
if and only if (II. 8p holds. 

In our second result, which is a Liouville-type theorem, we give a 
complete characterization of all finite Morse index solutions (whether 
radial or not, whether positive or not). 



Theorem 1.3. Let u be solution to ( TTjl ) with finite Morse index. 

• Ifp e (l,p c (n)), p ^Ps{n), then u = 0; 

• If p = Ps( n ), then u has finite energy i.e. 

(An) 2 = f \u\ p+1 < +00. 
If in addition u is stable, then in fact u = 0. 
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Remark 1.4. According to the preceding discussions, Theorem \1.3\ is 
sharp: on the one hand, in the critical case p = ps(n), Guo, Li and 
one of the authors [16] have constructed a large class of solutions to 
(11.11) with finite energy. Since in this case ( - p ~^ n = p + 1, by a result of 
Rozenbljum [29], such solutions have finite Morse index. On the other 
hand, forp > p c {n), all radial solutions are stable (see [TT] [20] ). 

Remark 1.5. The above theorem generalizes a similar result of Farina 
[13] for the classical Lane-Emden equation. 

Now consider (11. ip when Q is a smoothly bounded domain of MJ 1 and 
supplement it with Navier boundary conditions: 

f A 2 u = X(u + If in n 

\ u = Au = on dfl, [ } 

where A > is a parameter. It is well known that there exists a critical 
value A* > depending on p and fl such that 

• If A G (0, A*), (II. 9p has a minimal and classical solution u\, 
which is positive and stable; 

• If A = A*, a unique weak solution, called the extremal solution 
u\* exists for (P\*). It is given as the pointwise limit u\* = 
lim At u x ; 

• No weak solution of (II .9p exists whenever A > A*. 

An outstanding remaining problem is the regularity of the extremal 
solution U\* . An application of Theorem 11.31 and standard blow-up 
analysis gives 

Theorem 1.6. If n < n p (equivalently p < p c {n)), the extremal solu- 
tion Ux* is smooth. 

More generally, 

Theorem 1.7. Assume n < n p (equivalently p < p c {n)). 

• Let Q be a smoothly bounded domain and u G Hq(Q) fl H 2 (Q) 
be a solution of (I1.9P of finite Morse index k G N. Then, u 
is smooth and there exists a constant C > depending only on 
k, N, Q,p such that 

||«|U°°(n) < C. 

In particular, any classical solution satisfies the above inequal- 
ity. 

• Let Q be any open set and u G Hq(Q) fl H 2 (Q) be a solution of 
(II. ip . Then, there exists a constant C > depending only on 
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fc, N, Q : p such that for every i < 3, 

\V l u\ < Cdist (ar,9n)"^ _i a.e.m 0, 

The first part of the above theorem is again sharp since the singular 
solution defined by fll.Tf) . (jl.6p is stable whenever n > n p . For such 
dimensions, one can still try to estimate the size of the singular set of 
solutions. 

Definition 1.8. A point x belongs to the regular set of a function 
u G L} oc {Vl) if there exists a neighborhood B of x such that u G L°°(B). 
Otherwise, x belongs to S, the singular set of u. 

By definition, the regular set is an open set. By elliptic estimates 
applied to (11. ip . u is smooth in its regular set. Now, we state the 
interior partial regularity for U\* . 

Theorem 1.9. Let n > n p and let u\* be the extremal solution to 
(11.91) . Then the Hausdorff dimension of its singular set S is no more 
than n — n p . Moreover, when n = n p , S is a discrete set. 

We now list some known results. We start with the analogous second 
order equation 

Au+ {uf^u = 0, in R n . (1.10) 

As mentioned earlier, Farina completely classified finite Morse index 
solutions (positive or sign-changing) in his seminal paper [TJ]. His 
proof makes a delicate use of the classical Moser iteration method. 
More precisely, if one multiplies the equation (ll.lOp by a power of u, 
say u q , q > 1, Moser's iteration works because of the following simple 
identity 

/ u«(-Au) = j-^— 2 f |V^| 2 , V^GC^l-). 

There have been many attempts to generalize Moser's iteration tech- 
nique (or Farina's approach) to fourth order problems like ( 11. ip . Un- 
fortunately, this runs into problems: the corresponding identity reads 

/ n g (A 2 n) = 4g f \Au^\ 2 -q(q-l) 2 [ n 9 " 3 | Vn| 4 , Vn G C 2 (W ri ) 

and the additional term J Rn u q ~ 3 \ Vm| 4 makes the Moser iteration argu- 
ment difficult to use. 

Another strategy is to use the test function v = —Au. This allows 
to treat exponents less than + e n for some e n > 0, see the works 
of Cowan- Ghoussoub-Esposito [3] and Ye and one of the authors [35] . 
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Another approach, obtained by Cowan and GhoussoutQ[4j , and further 
exploited by Hajlaoui, Harrabi and Ye [TO], is to derive the following 
interesting interpolated version of the inequality: for stable solutions 
to (11. ip . there holds 



This approach improves the first upper bound + e n , but it again 
fails to catch the optimal exponent p c {n) (when n > 13). It should be 
remarked that by combining these two approaches one can show that 



stable positive solutions to (II. ip do not exist when n < 12 and p > - 



n-4' 



see [19] . 

In the above references, only positive solutions to (II. ip are consid- 
ered. One reason is their use of the following inequality, due to Souplet 

era 

2 i 

Au + m^ 1 < in 1" (1.11) 

As observed in [11] for a similar equation, the use of the above inequal- 
ity can be completely avoided. 

In this paper we take a completely new approach, which also avoids 
the use of (II. lip and requires minimal integrability. One of our motiva- 
tions is Fleming's proof of the Bernstein theorem for minimal surfaces 
in dimension 3. Fleming used a monotonicity formula for minimal 
surfaces together with a compactness result to blow down the mini- 
mal surface. It turns out that the blow-down limit is a minimal cone. 
This is because the monotonic quantity is constant only for minimizing 
cones. Then, he proved that minimizing cones are flat, which implies 
in turn the flatness of the original minimal surface. 

At last, let us sketch the proof of Theorem 11.31 we first derive a 
monotonicity formula for our equation (II. ip . Then, we classify sta- 
ble solutions: this is Theorem 14.11 in Section HI To do this, we esti- 
mate solutions in the L p+1 norm, utilizing the afore-mentioned meth- 
ods available in the litterature, and then show that the blow-down limit 

4 

u°°(x) = liniA^oo \ p ~ 1 u(\x) satisfies E{r) = const. Then, Theorem ll.il 
implies that u°° is a homogeneous stable solution, and we show in The- 
orem [3J] that such solutions are trivial if p < p c {n). Then similar to 
Fleming's proof, the triviality of the blow-down limit implies that the 
original entire solution is also trivial. Finally, we extend our result to 
solutions of finite Morse index. 



1 a similar method was first announced in [3] , and later published in the work by 
Farina-Sirakov and one of the authors |10) . 
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2. Proof of the Monotonicity formula 
In this section we derive a monotonicity formula for functions u G 



W 4 ' 2 {B R {0)) n L p+l {B R {0)) solving flUD in B R (0) C Q. We assume 
that v>— A - 

r n—A 



Proof of Theorem \l.l[ Since the boundary integrals in E(r;x,u) only 
involve second order derivatives of u, the boundary integrals in ^p(r; x, u) 
only involve third order derivatives of u. By our assumption u G 



W^ 2 {B R {0))nL p+1 (B R (0)), for each B r (x) C B R (0), u G W 3 ' 2 (dB r (x)). 



Thus, the following calculations can be rigorously verified. Assume that 
x = and that the balls B\ are all centered at 0. Take 



u x (x) := \p- 1 u(\x), v x (x) := \p- 1+2 v(\x). 
We still have v x = Au x , Av x = (u x ) p , and by differentiating in A, 

. du x dv x 

A = — . 

dX dX 

Note that differentiation in A commutes with differentiation and inte- 
gration in x. A rescaling shows 




Define 



v := Am 



and 




Hence 



dX 




(2.1) 
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x d du x dv x du x 
v x 



<dBi 9r dX dr dX 

In what follows, we express all derivatives of u x in the r = |x| variable 
in terms of derivatives in the A variable. In the definition of u x and v x , 
directly differentiating in A gives 

du x , . 1/4 X/ . du x , A , . 

luW = A(^T" W + r ^ W )' (2 - 2) 

dv x , , 1 /2(p + l) A/ <9u A , A . , 

*rW = 5(T=T V(l)+r -* (l) )- (2 - 3) 

In ( 12.21) . taking derivatives in A once again, we get 

,d 2 u x , , cfot A , , 4 d« A , d du x , . , 

x H>F {x) + lx ix) = —iHx {x) + r TrHx {x) - {2A) 

Substituting (Q and fl2l| into (JE!) we obtain 

d£ /" ^ x fd 2 u x + p - 5 du x \ du x fdv x 2(p + l^ 

^ 7aBi V ^ 2 p — I dX J dX \ dX p 

Xv \ d 2 u x 2v x — - X 
ob 1 dX 2 dX 

Observe that v x is expressed as a combination of x derivatives of u x . 
So we also transform v x into A derivatives of u x . By taking derivatives 
in r in ( 12. 2 p and noting ( 12. 4p . we get on dBi, 

d 2 u x d du x p + 3 du x 

X- 





du x 


)- 


dX 


du x 


dv x 


dX 


dX ' 



dr 2 dr dX p — 1 dr 

„d 2 u x p-5,du x p + 3fdu x 4 , 
= X 2 —— + - — -A— — - ( A— -u x 



dX 2 p — 1 dX p — 1 \ dX p — 1 
2 d 2 u x 8 du x 4(p + 3) A 

~dxJ'p~^i ~dx + {p-iy u ' 



Then on 8B\ 



, d 2 u x n — 1 <9w A 1 . x 
u = -<rr + ^ + ^ A ^ 

, 2 d 2 u x 8 du x 4(p + 3) \ , ^ ( ,du x 4 A 
x2 d 2 u x ( 8 \ du x 4,4 , , A , 
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Here Ag is the Beltrami-Laplace operator on dBi and below repre- 
sents the tangential derivative on dB\. For notational convenience, we 
also define the constants 

8 4 4 

a = n-l -, (3 = --n+2. 

p — 1 p — 1 p — 1 

Now reads 

d /• . ( x2 d 2 u x ^du x A d 2 u A 

+3(A 2 ^ + aA^ + /^ dtt * 



c/A 2 c/A / dA 

, du x d ( , 9 d 2 u A , du A „ , 

/" \ d 2 « A . , d« A , d« A . d« A 

+ / AA^_ + 3 A^— -A— A 9 — 



dB 



— R± + R 2 . 

Integrating by parts on <9i?i, we get 

/" x „ d 2 u A A „ du x , ,„ dw A , 2 

dA 2 dA 1 dA 

X d 2 f f ,„ A2 \ 3 d f f A A /" du 



2 dA 2 / 2 dA / J,^ dX 

For after some simplifications we obtain 

f , A 2 ^ A ,dw A „ A d 2 w A 
ifc = / A A 2 — — + aX— + f3u } 



dB! 



dX 2 dX ' ) dX 2 

„d 2 u x , dw A „ A d« A 
A 2 —— + aA— + f3u x 



dX 2 dX J dX 

du x f. 2 d 3 u x .d 2 u x . .du x 



dA 2 / rfA 2 tiA dA 2 ^ dA 
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+(2a-B)X 
' d 2 u 



du x 



2X d 



dB 1 



B d 2 



dX 2 
AM 



dX t 
+ 4A 2 



, du x d 3 



dX dA 3 
d 2 u x du x 



dX 2 dX 



+ (2a - 2/3) A 



2dA 2 

Here we have used the relations 



1 d 
2dX 



d ( du 



an 2 



dX V dA 



dw A 

/3 d 
2 dA 



(« A ) : 



A 



/ 2 -2ff-X(f) 



'\2 



and 



-A 3 /T 



A 3 



((f) 2 )' 



+ 3AYr + A 3 (f) 2 



Since p > ^4, direct calculations show that 



n-4> 



a — /3 
Thus, 



n 



2X 3 



p — 1 



4 



4 



2A A 



dX 2 
d 2 u x 



+ 4A 



p — 1 \p — 1 

,d 2 w A d« A 



n + 2 > 1 



dA 2 dA 



+ (2a - 2/3) A 



du 



AN 2 



dA 2 dA 



+ (2a - 2/3 - 2)A 



du 
~dA~ 



du 
~d~X 

AN 2 



AN 2 



(2.6) 



(2.7) 



> 0. 



Then, 
Ri > 



B d 2 
2d>? 



A^ 2 



A («») 



1 d 
2dA 



d / du A A 2 



dA V dA 



B d , XN 2 



Now, rescaling back, we can write those A derivatives in R\ and R2 as 
follows. 



d 2 



d ( A\2 d / 8 + i_ n /" 5 

5a (« ) = 5a ( A '-' L" 



asi ^A 2 



A\ 2 



A (« ) 



dA 2 



A?" 1 



9Bi 



dA 



d f du 



.AN 2 ' 



dA 2 



A 



dA V dA 

A|2 



d 

dA 



d 



dA 



9Bi 



dA 2 



A l+-| T+ 2 + l-n 



4 



(9m , 



9B\ 
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dX 



~TT I / \V 6 U X \ 2 I = — 



A ^T+2+l-n 



9B A 



Vm| 2 - I — I 
or 



Substituting these into -4rE(X; 0, u) we finish the proof. 



□ 



Denote c(n,p) = 2a - 2(3 - 2 > 0. By fl2"77jl . we have 
Corollary 2.1. 

rf by n ^ ^ /- ^ -"+2+^ f ( 4 _ x <9u 
— -fr(r; U, w) > c[n,p)r p- 1 / r w + — 

dr J 9Br VP - 1 

In particular, if E(X; 0, w) = const, for all A G (r, i?) , u is homogeneous 
in -Br \ -B r : 



9« v 



M (x) = \X\ p- x U 



\x\ 



We end this section with the following observation : in the above 
computations we just need the inequality (12.61) to hold. In particu- 
lar the formula can be easily extended to biharmonic equations with 
negative exponents. We state the following monotonicity formula for 
solutions of 



A u = , u > in Q C R n . 

vP 



Lemma 2.2. Assume that p satisfies 
n-2 



4 4 
— r( + n-2) 



(2i 



(2.9) 



p+1 p + I p + 1 

Let u be a classical solution to Ii2.8\) in B r (x) C Br(x) C Q. Then the 
following quantity 



E(r: x, u) 



f p+i 



1 



(An) 



B r (x) 



P — 1 



-U 



1-p 



— — (n-2-\ — )r p+ 1+1 n I 

P+1 V P+ 1/ J9B r (x) 

n — 2 H I — I r P+ 1 / u 

dB r (x) 

du\' 



p+1 

r 3 d 
+ lldr 



p + 1 ) dr 



f p+i 



-l-n 



9B r (x) V P+1 



r u + — — 



<9r / 



+ 



1 d 
2~d~r 



p+i 



dB r (x) 



v« 2 - 1— | 2 
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2 J d B r (x) V 9r J 

is increasing in r. Furthermore there exists Co > such that 

^M,«)> w —- - /^(--l r -H, + |) (2,0) 

In the rest of the paper, sometimes we use E(r;x) or E'(r) if no 
confusion occurs. 

3. Homogeneous solutions 

For the applications below, we give a non-existence result for homo- 
geneous stable solution of (11. ip . (This corresponds to the tangent cone 
analysis of Fleming.) By the Hardy- Rellich inequality, this result is 
sharp. 

Theorem 3.1. Let u G W^(K n \ {0}) be a homogeneous, stable 
solution of (II. ip in M. n \ {0}, for p G (^zf ,p c (^))- Assume that 
\u\ p+1 G L} oc {R n \ {0}). Then u = 0. 

Proof. There exists aw 6 W 2,2 (S n ~ 1 ) such that in polar coordinates 

4 

w(r, 9) = r p- 1 w(9). 

Since u G W 2 ' 2 (B 2 \B 1 )nL p+1 (B 2 \B 1 ), w G ^(S™" 1 ) nL p+1 (S n ~ 1 ). 
Direct calculations show that w satisfies (in W 2,2 (E> n ~ 1 ) sense) 

A 2 9 W - JxAgW + J 2 W = w p , (3.1) 

where 

J 1= (-*- + 1 ) („-4— Q+ — f„-2— L. V 

Vp-i y v p-iy p-i\ p-iy 
j2 p - 1 (p - 1 + 2 ) ( n 4 p - 1 ) ( n 2 p - 1 ) ' 

Because u> G iy 2 ' 2 (S n ~ 1 ), we can test ( 13. ip with w, and we get 

/ \A e w\ 2 + J 1 \V e w\ 2 + J 2 w 2 = [ \w\ p+1 . (3.2) 

For any £ > 0, choose an r/ e G C^°((|, |)), such that r] e = 1 in (e, i), 
and 

r l ? 7e( r )l + r ' 2 | ? 7e( r )l — 64 for all r > 0. 

Because w G H /2 ' 2 (S™~ 1 ) nL p+1 (§ n_1 ), r~ !1 ^ ± w(6)r] £ (r) can be approxi- 
mated by CH5 4/£ \ B e/4 ) functions in W 2 ' 2 ( J B 2/e \ S e/2 ) n U+\B 2/e \ 
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B e /2). Hence in the stability condition for u we are allowed to choose 
a test function of the form r~ R ^w(9)rj £ (r). Note that 

A (j-~ !L ^ ± w(6)r) E (r) S j = — n ^ n ^ — r~%r) e (r)w(9) + r~%r] £ (r)A e w(9) 

_3 r -f + 1 rj ' e ( r ) w (e) + r~^ +2 ri"(r)w(6). 
Substituting this into the stability condition for u, we get 

p \w\ p+1 d9^J f jf r-\(r) 2 dr 

/ /-+00 N 

+0[l / rr/^(r) 2 + r 3 ?7^(r) 2 + r 2 |?7g(r)|r7 £ (r) + r?7 e (r)|?7^(r)|(ir 

x ( I w(9) 2 + \V e w(9)\ 2 de)}. 
Vis™- 1 / 

Note that 



p+oo 

I r~ lr q £ (r) 2 dr > | log erj 
Jo 



+00 

rr]' e (r) 2 + r 3r q"(r) 2 + r 2 \r]' £ (r)\r] £ (r) + rrj e {r)\r}"{r)\dr < C, 
for some constant C independent of e. By letting e — > 0, we obtain 

p [ \wr- x dB< [ |A^| 2 + ^fil|V^p , 

Jgn-l J§"-1 2 tO 

Substituting ( 13. 2 p into this we get 

/ (p-i)|A^| 2 +(M-^^)|V^| 2 +(M-^^V 2 < 0. 

If=±£ <p<p c (n),thenp-l >0,pJ 1 -^^ > and pj 2 - " 2 ^~ 4 ) 2 > 
(cf. p. 338 in [T8]), so iy = and then u = 0. □ 

For applications in Section 6, we record the form of E(R; 0,u) for a 
homogeneous solution u. 

4 

Remark 3.2. Suppose u(r,0) = r p~ 1 w(9) is a homogeneous solution, 
where p > £±| and w G W 2,2 ^™" 1 ) D L p+1 (S n_1 ). in to case, /or any 
r > 0, 

I A 1 2 1 I |p+l ^ n-4^±l 

I Am I + |m| < cr p- 1 . 



r2 2 (n-4) 2 9 
-to 2 . 



B r \B r / 2 
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Because n — 4^j- < 0, by choosing r = 2~ l R and summing in i from 1 
to +oo, we see 



Ib, 



\ Au \ 2 + \u\ p+1 < cR n ~^ . 



'B R 

which converges to as R — > 0. Hence for any R > 0, E(R;0,u) is 
well-defined and by the homogeneity, it equals E(l;0,u). By definition 

E(l;0,u) = [ \{Auf - -^-\u\ p+l 

Bt * P+1 



4 / 4 



P~ 1 V P 1 / JdB 1 JdB 1 

1 1 \ /" . 1 /" f du . dAu\ 



\u\ 



p+1 - "" Au _ u 



2 P+1/ 7bi 2 VaBi V^r / 

n-2 i-] / m 2 + / |V e w| 2 . 



p — 1 \ p — 1 

5?/ noting that 

du 4 _ 1 <9 2 u 



SBi JdB! 



r L u, — — - = I + 1 ) r 2 u 



dr p — 1 dr 2 p — 1 \p — 1 



— — = - 2 H r -1 Au, Am = h 2 - n r~ A u+r~ 2 A e u, 

ar \ p — 1 / p — 1 \p — 1 / 

we get 

Replacing \u\ p+1 by (Au) 2 , we also have 

+— i- (n-2 ^- J / it 2 . 

P + 1 V P- 1/ 73B! 

4. The blow down analysis 

In this section we use the blow-down analysis to prove the Liouville 
theorem for stable solutions. Throughout this section u always denotes 
a smooth stable solution of (II. ip in W 1 . 

Theorem 4.1. Let u be a smooth stable solution of (11. ip on W 1 . If 

1 < p < p c (n), then u = 0. 

The following lemma appears in [35] for positive solution. It remains 
valid for sign-changing solutions, see also [19] . 
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Lemma 4.2. Let u be a smooth stable solution of (II. ip and letv = Au. 
Then for some C we have 

[ (v 2 + M p+ V <C [ u 2 (|V(At7) ■ V77I + (At/) 2 + |A(|Vr?| 2 )|) dx 

JR n JR n 



+C j |m;||V?7| 2 (ix 

JR n 

(4.1) 

for all rj £ C£°(R n ). 

Proof. For completeness we give the proof. We have the identity 

(A 2 ^V 2 dx= [ (A(^)) 2 + / (-4(Ve-Vr / ) 2 + 2^|Vr/| 2 )^ 



+ f e 2 (2V(Ar/) • V77 + (A V ) 2 )dx, 

JR™ 

for £ £ C 4 (lR n ) and 77 £ Co°(M n ), see for example Lemma 2.3 in [35] . 
Taking £ = it yields 

/ |u| p+ V^= / (AN) 2 + / (-4(Vm- Vr/) 2 + 2wt;|Vr7| 2 )rfx 
Jr™ Jr™ Jr™ 

+ / u 2 (2V(Ar]) ■ Vr] + (Ari) 2 )dx, 

JR n 

Using the stability inequality with ur\ yields 

V I \u\ p+1 r) 2 dx < j (A(un)) 2 . 

JR n JR n 

Therefore 

/ (\u\ p+1 r} 2 + (A(ur])) 2 )dx < C [ (\Vu\ 2 \Vri\ 2 + \uv\\Vri\ 2 )dx 

JR n JR n 

+ C I u 2 {\ V(Ar}) ■ Vr]\ + (Arj) 2 )dx. 

JR n 

Using A(r)u) = vrj + 2Vr] ■ Vu + 77 Am we obtain 

(\u\ p+1 + v 2 )r} 2 dx < C [ (|Vw| 2 |Vt7| 2 + \uv\\Vr]\ 2 )dx 

l n JR n 

+ C f u 2 (\V(Ar])-Vr]\ + {Ar]) 2 )dx. 

JR n 



But 



2/ \Vu\ 2 \V7]\ 2 dx= I A(u 2 )\Vr]\ 2 dx -2 I uv\Vr]\ 2 dx 

JR n JR n JR n 
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= / u 2 A(\W7]\ 2 )dx - 2 / uv\X7r]\ 2 dx, 

and hence 

(M p+1 + v 2 ) V 2 dx < C [ u 2 (|V(Ar/) ■ V77I + (Ar/) 2 + | A(|Vr/| 2 )|)dx 

+ C7 ! \uv\\Vri\ 2 )dx. 

This proves flUD □ 
Corollary 4.3. There exists a constant C such that 

[ v 2 + \u\ p+1 < CR- 4 [ u 2 + CR- 2 [ \uv\, 

JB R (x) JB 2R (x)\B R (x) J B 2R {x)\B R (x) 

(4.2) 

and 

v 2 +\u\ p+1 <CR n ~ 4 ^. (4.3) 

B R (x) 



for all Br{x). 



Proof. The first inequality is a direct consequence of ( 14.11) , by choosing 
a cut-off function r\ G (B 2 r(x)) , such that 77 = 1 in B R (x), and for 
A; < 3, \V k r]\ < 

Exactly the same argument as in [35] or [19] provides the second 
estimate. For completeness, we record the proof here. Replacing 77 in 
(14.11) by 77 m , where m is a larger integer and 77 is a cut-off function as 
before. Then 

| m ;||Vr7 m | 2 = m 2 [ Mr/ 2m - 2 |V7?| 2 

J B 2R (x)\B R (x) 
1 ' 22m 1 ri f „.2„2m-4|v7„|4 



< ^ / vY m + C I «V |Vr?| 



Substituting this into (14. ip . we obtain 

y + | u |f+ 1 ) 7? 2 - < cir 4 / n 2 r/ 2m - 4 



< cr- 4 [ / | M ry m - 2 ^ +i M 

'-Bail W 



2 

p+1 



_2_ 

P+1 



< C7T 4 / | u |P+ym-2)CP+l) ^-PTT). 

\Jb 2R {x) J 

This gives ( I4.3p . Here we have used the fact 77 2m > r^™ -2 )^ 1 ) because 
< i] < 1 and m is large. □ 
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Proof of Theorem \4.1\ for 1 < p < For p < we can let R — > 
+00 in (14.31) to get u = directly. If p = this gives 



So 



v 2 + \u\ p+1 < +00. 



lim / v 2 + M p+1 = 0. 



Then by (14. 2p . and noting that now n = 4|±j, 



„2 + < CjR -4 / U 2 + C \V\ 2 

Br(x) JB 2r (x)\B r (x) Jb 2 r(x)\B r (x) 



2 

< C7T 4 ( / \u\ p+l Y +1 R n(l -^ + C 



\v\ 



B 2R (x)\B R (x) / J B 2R {x)\B R {x) 



< C / +c 

\J B 2R (x)\B R (x) J JB 2R (x)\B R (x) 

This goes to as R — > +00, and we still get u = 0. □ 
Next we concentrate on the case p > We first use (14. 3 p to show 
Lemma 4.4. lim E(r; 0, u) < +00. 

r—t+oo 

Proof. Since E(r) is non- decreasing in r, we have 

1 p2r 1 p2r rt+r 

E{r) <- E(t)dt < — / / E(X)d\dt. 

r Jr r J r Jt 

By <m, 




r Jt 



r 

Next 

^ f2r r-t+r / g 




AA +1 '" / m 2 I dXdt 

r Jt JdB x 



1 /"2r 




n u(x) 2 dxdt 



r JB t+r \B t 



< - 2 f r {l M^ +1 -^V +1 ( [ \u(x)r^~a 

r Jr \JB 3r \ Br J \JB 3r 

< C. 
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The same estimate holds for the term in E(r) containing 



dB x 



v« - bH 

or 



For this we need to note the following estimate 



V< < Cr 2 I (Auy+Cr 

B 2r 



13 , 

Now consider 

2r r t+r A 3 rf 



-2+n 



£-1 
P+l 



p+l 



B 2 r 



2 
P+l 



1 

2r^ 



2r 



3 

2^2 



Ap-i 

2 dX 

{(t + r)^ +4 ~ n 

2r />t+r 



-1-n 



4 fl»V 
p — 1 <9r y 



dXdt 



4 . . i duY 

2 



4 



I*"'" ' <9r ) 

dB t \P — 1 ar J 

4 



r Jt 



< 



c 



< c. 



\x 



\ -2-7+3-71 

-^r+4-n 

p-l 



9B\ 



}dt 



du\' 



X 1 u + —) dXdt 
p — 1 or J 



Bsr\Br 



P — 1 



i ,_i du\' 

\ x \ u + ^~ 
ar / 



□ 



The remaining terms in E(r) can be treated similarly. 
For any A > 0, define 

u x (x) := A~n(Ax), v x (x) := X~ +2 v(Xx). 

u x is also a smooth stable solution of (II. ip on R n . 

By rescaling (TOD, for all A > and balls B r .(x) C R n , 

( w A ) 2 + | n A| P+ l< C , r n-4^i_ 

In particular, u x are uniformly bounded in Lf |j 1 (]R n ). By elliptic es- 
timates, u x are also uniformly bounded in Wf^iW 1 ). Hence, up to a 
subsequence of A — > +oo, we can assume that u x — > u°° weakly in 
Wf£ (M n ) D L^^M"). By compactness embedding for Sobolev func- 
tions, u x — > u°° strongly in W Zo ' c (M n ). Then for any ball B R (0), by 
interpolation between L q spaces and noting (14.31) . for any q G 
as A +oo, 



A „ oo i 



I«"-«~IMb h (o)) < ll^-^llli^flW)!!^-^!!^*!^^)) ( 4 - 4 ) 



A . oo ii 1— t 
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where t G (0, 1) satisfies ± = t + ±=f . That is, u A -»■ u°° in Lf oc (M n ) for 
any g G + 1). 

For any function y? G (M n ), 

Au°°Av? - (u 00 )^ = lim / Au A A<^ - (u A )V = 0. 

(Av9) 2 -p(n°°) p "V 2 = lim f (Ay?) 2 — p(u x ) p ~ 1 ip 2 > 0. 

Thus u°° G M^ 2 (E n ) n ifi^R") is a stable solution of (JTTJ) in R n . 
Lemma 4.5. u°° is homogeneous. 

Proof. For any < r < R < +oo, by the monotonicity of E(r; 0, it) 
and Lemma [4.41 

lim E{XR;0,u) - E(\r;0,u) = 0. 

A— >+oo 

Then by the scaling invariance of E and applying Corollary 12.11 we see 
= lim E{R; 0, u x ) - E(r; 0, u x ) 

A— >+oo 



^|x|-V(x) + ^(x)) 



2 



> c(n,p) lim / — g — dx 



*->+°o .1 b r \b t \x\ n 2 ~ p- 1 



2 



> c(n,p) / ^_ dx. 

Note that in the last inequality we only used the weak convergence of 
u x to u°° in Wi^(R n ). Now 

4 _ 1 „ 9m 



p — 1 dr 
Integrating in r shows that 



r~ l u°° + — — = 0,a.e. in M n . 



u (z) = \x\ p-^u 



That is, m°° is homogeneous. □ 

By Theorem 13.11 u°° = 0. Since this holds for the limit of any 
sequence A — > +oo, by (14.41) we get 



lim u x = strongly in L 2 (B 4 (0)). 

A— »+oo 



Now we show 



Lemma 4.6. lim E(r;0,u) = 0. 

r— »+oo 
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Proof. For all A — > +00, 

lim / [u x f = 0. 

Because t> A are uniformly bounded in L 2 (I?4(0)), by the Cauchy in- 
equality we also have 



lim / |uV| < lim ( / (u x ) 2 



v x ) 2 



By (Q, 

lim / (v x ) 2 + \u x \ p+1 < C lim 



(u x ) 2 + / \u x v x 



A ^+°°JB 3 (0) A^+oo \Jb 4 (o) -/b 4 (o) 

= 0. (4.5) 

By the interior L 2 estimate, we get 



lim 

A— >+oo 



/ V|W| 2 = o. 



'^(0) k<2 

In particular, we can choose a sequence Aj — > +00 such that 



/ £|VVf<2" 



By this choice we have 

„2 +00 „ +00 „2 



£ / £ I Wf rfr < £ / / £ I Wf rfr 

i=l "^-Br fc<2 i=l "' 1 "' l9 ' B '- fc<2 



< 1. 



That is, the function 

+00 „ 

/ ElWf €^((1,2)). 

j=l fc<2 

There exists an r G (1, 2) such that /(r ) < +00. From this we get 

lim \\u Xi \\ W 2,2( dBrn ) = 0. 

Combining this with (14.51) and the scaling invariance of E(r), we get 
lim E{Xir ;0,u) = lim E(r Q ; 0, u Xl ) = 0. 

i— >+oo i— >+oo 

Since Ajr — >■ +00 and £'(r; 0, u) is non- decreasing in r, we get 

lim £(r;0,w) = 0. □ 

A— >+oo 
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By the smoothness of u, ]xmE(r;0, u) = 0. Then again by the 

r-i-0 

monotonicity of E(r; 0,w) and the previous lemma, we obtain 

E(r; 0, u) = for all r > 0. 

Then again by Corollary 12. 1[ u is homogeneous, and then u = by 
Theorem 13.11 (or by the smoothness of u). This finishes the proof of 
Theorem 14.11 

5. Finite Morse Index Solutions 

In this section we prove Theorem 11.31 First, by the doubling lemma 
[25] and our Liouville theorem for stable solutions Theorem 14.11 we 
have 

Lemma 5.1. Let u be a finite Morse index (positive or sign changing) 
solution of ( Eip . There exists a constant C\ and R such that for all 

x e B Ro (o) c , 

4 

< C\x\ p- 1 . 

Proof. Assume that u is stable outside Br q . For x G B c Ro , let M(x) = 

\u(x)\~ and d(x) = \x\ — R , the distance to Bj^. Assume that there 
exists a sequence of Xk 6 such that 

M(x k )d(x k ) > 2k. (5.1) 

Since u is bounded on any compact set of M™, d(x k ) — > +oo. 

By the doubling lemma [25], there exists another sequence y k G Bj^, 
such that 

(1) M(y k )d(y k ) > 2k; 

(2) M(y k ) > M(x fc ); 

(3) M(z) < 2M{y k ) for any z G B c Rq such that \z - y k \ < j^- y 
Now define 

u k (x) = M(y k )~^u(y k + M^^x), for x G B k (0). 

By definition, |u fc (0)| = 1. By (3), \u k \ < 2^ in B k (0). By (1), 
B k /M{ yk )(yk) n B Ro = 0, which implies that u is stable in B k/Miyk) (y k ). 
Hence u k is stable in -Bfc(O). 

By elliptic regularity, u k are uniformly bounded in C 5 (B k (0)). Up to 
a subsequence, u k converges to in C^ c (lR n ). By the above conditions 
on u k , we have 

(1) |uoo(0)| = l; 

(2) \ Uoo \ < 2V in W 1 ; 

(3) is a smooth stable solution of fll.ip in M n . 
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By the Liouville theorem for stable solutions, Theorem 14.11 Uoo = 0. 
This is a contradiction, so ( 15.1 p does not hold. □ 

Corollary 5.2. There exists a constant C\ and R such that for all 
xeB 3Ro (0) c , 

^|a:|^ +fc |V fc n(x)| < C7 3 . (5.2) 

fc<3 

Proof. For any Xo with \xq\ > 3Ro, take A = ™ and define 

4 

u{x) = \p- 1 u(x + Ax). 

By the previous lemma, \u\ < C\ in -Bi(O). Standard elliptic estimates 
give 



5^|V fc n(0)| <C7 3 



fc<3 



Rescaling back we get (15. 2p . □ 

Remark 5.3. By the same proof of Lemma \5. 1] and Corollary \5.^ one 
easily obtains the second part of Theorem \l.l\ 



5.1. The subcritical case 1 < p < We use the following Po- 



rt— 4 



hozaev identity. For its proof, see [26 
Lemma 5.4. 

^(A^-^-M^ (5.3) 
b r 2 p+1 

R,. , 9 R . ^dudAu . d(x-Vu) 

-(An) 2 + — -\u\? +1 + Rj-^— - Au-±- >-. 

9Br 2 p+1 or or or 

By taking R — > +oo and using (15. 2p . and noting that p < we 
see that 

R.. , 2 R . .„,-, dudAu . d(x-Vu) 

By (15. 2p . we also have 

(Au) 2 + \u\ p+1 <C7(l + |x|)- 4 ^. 
Since p < ^ , 4^ > n. Hence 

^ n— 4' p— 1 



/ (Am) 2 + M p+1 < +oo. 



Taking limit in ( I5.3p . we get 



71 — 4 ... „ // 
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Take an 77 G C£°(-B 2 ), 77 = 1 in £>! and J^ fc<2 |V fc r?| < 1000, and denote 
Vr( x ) = ^{x/R). By testing the equation ( 11. ip with u(x)ri 2 R) we get 

/ (Au) 2 V 2 r-\u\ p+ W r = - [ 2VuV V 2 R + uArj 2 R . 

By the same reasoning as above, we get 

/ (Am) 2 - \u\ p+1 = 0. 

Substituting (15.41) into this, we get 
n — 4 ri 



2 p 



= 0. 



Since - ^ < 0, u = 0. 

5.2. The critical case. Since ■» is stable outside Bf> . Lemma 14.21 still 
holds if the support of 77 is outside B Ro . Take <p G C^°(B 2 b. \ B2Rq), 
such that </? = 1 in £> R \ -B3_r and ^ fc<3 |x| fc |V fe (^| < 100. Then by 
choosing 77 = ip m , where m is large, in (14. ip . and by the same reasoning 
to derive ( 14. 3p . we get 



/ 

Jb, 



(Atj) 2 + \u\ p+1 < C. 

'B R \B 3R() 

Letting R — > +00, we get 

/ (Aw) 2 + \u\ p+1 < +OO. 

This then implies that 

lim / R- 1 \Vu\ + R- 2 \u\ = 0. 

Then we can proceed as in the subcritical case to prove that 

! (Am) 2 - \u\ p+1 = 0. 

5.3. The supercritical case. Now we consider the case p > 

Lemma 5.5. There exists a constant C 2 , such that for all r > 3Rq, 
E(r;0,u) <C 2 . 

Proof. Expanding those boundary integrals in E(r; 0, u) into a full for- 
mulation involving the differentials of u up to third order, and substi- 
tuting (15. 2p into this formulation, we get 



(7 (A«) 2 4 


~ I«| P+1 ) 


\J B r 





p+1 ) + Cr^ +1 - n [ u 2 

J JdB r 
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+Cr^ +2 ~ n [ \u\\Vu\ + Cr^ +3 ~ n [ \Vu\ 2 

JdB r JdB r 

+CV^ +4 " n [ |Vu||V 2 u| 

JdB r 

< c. 

This constant only depends on the constant in ( 15. 2p . □ 

By Corollary l2.lt we get 
Corollary 5.6. 



n 1 8 

B c \r\ p-i 



dx < +oo. 



As in the proof for stable solutions, define the blowing down sequence 

, 4 

u [x) = \p- 1 u(Xx). 



By Lemma 15.1 1 u are uniformly bounded in C (B r (0) \ Bi/ r (Q)) for 
any fixed r > 1. u x is stable outside -Br o /a(0). There exists a function 
u°° G C 4 (lR n \ {0}), such that up to a subsequence of A — > +oo, u x 
converges to u°° in C^ c (R n \ {0}). u°° is a stable solution of ( II. ip in 
R n \ {0}. 

For any r > 1 , by Corollary 15. 6| 

2 



m-2- 



B r \B 1/r \x\ p 

^ixrv(x) + ^(x)) 2 

lim / — o — dx 



A^+oo .l Br \B 1/r \ X \ n 2 ~ p- 1 

^-r\x\- l u(x) + " 

lim / — 5 — dx 



A^+oo J Bxr \B x/r \ X \ n 2 ' P-i 

= 0. 

Hence w°° is homogeneous, and by Theorem 13.11 u°° = 0. This holds 
for every limit of u x as A — > +oo, thus we have 

4 

lim |x| p- 1 u(x)| = 0. 

x— >oo 

Then as in the proof of Corollary 15.21 we get 
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Lemma 5.7. 

lim V|x|^ +fc |V fc w(x)| = 0. 

x— >oo ' ' 

For any e > 0, take an i? such that for \x\ > R, 
^|x|^ +fc |V fc n(a;)| < e. 

fc<4 

Then for r>i?, 

E(r;0,u) < Cr 4 ^- n [ [(Au) 2 + \u\ p+1 ] + Cer 4 ^~ n [ 

JB R (0) JB r (0)\B R (0) 
n 4 £±i +1 _ n f |-4£±i 

+Cer p- 1 / \x\ p- 1 

JdB r (0) 

< C(R)(r^~ n + e). 

Since A^^—n < and e can be arbitrarily small, we get lim E(r; 0, u) = 
0. Because lim E(r; 0, u) — (by the smoothness of u), the same argu- 

r— >0 

ment for stable solutions implies that u = 0. 

Remark 5.8. T/ie monotonicity formula approach here is in some 
sense equivalent to the Pohohazev identity method (see for example 
[35] ). The convergence of u x can also be seen by writing the equation 
in exponential polar coordinates. 



6. Partial regularity in high dimensions 

Here we study the partial regularity for the extremal solution to the 
problem (jl.9p , and prove Theorems 11.61 and 11.91 Recall that we defined 
n p to be the smallest dimension such that Theorem 13.11 does not hold. 
This is also the smallest dimension such that the Liouville theorem for 
stable solutions, Theorem 14.11 and the classification result for stable 
homogeneous solutions, Theorem 13 .![ do not hold. 

Proof of Theorem \l.b\ For < A < A* let u\ > be the minimal 
solution of ( 11. 9ft . We claim that 

SUp ||MA|U°°(fi) < +00. (6.1) 
Ae(0,A*) 

Then by elliptic estimates, as A — > A*, U\ are uniformly bounded in 
C 5 (f2). Because u\ converges to u\* pointwisely in Q, u x * £ C 4 (f2), 
and then we get u\* £ C°°(fl) by bootstrapping elliptic estimates. 
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To prove (16.11) . we use the classical blow up method of Gidas and 
Spruck. Let x x attain max.^u x , and assume that 

L\ = u x (x x ) + 1 +oo. 
By the maximum principle, x x G Q is an interior point and 

- Au x > in Q. (6.2) 

Define 

u x = L^ 1 (u x (x x + L x ~x) + lj in Q x , 

where £l x = L x 4 (Q — x x ). u x is a smooth stable solution of (II. ip in 
Q x , satisfying 

u x (0) = maxu x = 1, (6.3) 

and the boundary condition 

u x = L\~ \ Au x = on dfl x . 

From this, with the help of standard elliptic estimates, we see for any 
R > 0, u x are uniformly bounded in C 5 (Q X n -Br(O)). By rescaling 

- Au x > in Q x . (6.4) 

Since Q is a smooth domain, as A — > A*, Q x either converges to 
IR™ or to a half space H . In the former case, u x converges (up to a 
subsequence) to a limit u in C^ c (M n ). Here u is a positive, stable, C7 4 
solution of (II. ip in M n . Then by Theorem 14.11 u = 0. However, by 
passing to the limit in (16. 3p . we obtain 

«a(0) = I- 

This is a contradiction. 

If 17a converges to a half space H = {x\ > —h} for some h > 0, u x 
converges (up to a subsequence) to a limit m in Cf oc (H). Here m is a 
positive, stable, C 4 solution of (11.11) in H, with the boundary conditions 

u = Au = on OH. 

By taking limits in (16. 3p and (16. 4p . we obtain 

— Au — v > 0, in iJ, 

- Aw = n p > 0, in if, 
m(0) = maxii = 1. 

H 

By elliptic estimates, the last condition implies that v is bounded in 
H. Then by [Theorem 2, [6]] or [Theorem 10, [31]], Jg- > 0, g- > 0. 
Then the function w(y) = lim :) ; 1 _ >+OC) u(xi,y) exists for all y G IR™^ 1 and 
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satisfies A 2 w = w p in M™ -1 . By the arguments in Section 3 of [32] this 
function w must be stable in M. n ~ l and non trivial. By Theorem 11.31 
p > p c [n — 1) > p c (n). This is impossible. 

We conclude that u = 0, which is a contradiction. This finishes the 
proof of (JET). □ 



The remaining part is devoted to the proof of Theorem II .91 First we 
need the following lemma. 

Lemma 6.1. There exists a constant C , such that, for any ball B 2r {x) C 
r^~ n [ {u x * + l) 2p <Cr 4 ^ n [ (« A , + 1) P+1 + (A« A ,) 2 . (6.5) 



JB r {x) JB 2r {x) 

Proof. Denote w\ — U\ + 1. By the maximum principle and Lemma 
3.2 in [3], for any A G (0, A*), 



2A s+1 

Awj, < -4/ w x 2 < in Q. 

p + 1 

Since w\ is smooth in Q, we can follow the proof in [32] to get Eq. 
(2.15) in [35]. That is, for any rj G C °°(fi), 

[wf V 2 < [ -Aw x w p x {\V V \ 2 + \A V 2 \) (6.6) 

+C / (Aw A ) 2 [IVA77V77I + |A|Vr/| 2 | + |Ar/| 2 .] 
Jn 

Take </? G C^°(i?2r-(^)) such that < <p < 1, <£> = 1 in B r (x) and 

^r fc |VVl < 1000. 

fc<4 

Substituting 77 = ip m into (16.61) with m larger, and then using Holder 
inequality (exactly as in the derivation of Eq. (2.16) of [35]), we get 
(16.51) for u\. 

This implies that u\ are uniformly bounded in L 2p c (Q). By the in- 
terior I? estimate, u\ are also uniformly bounded in W^ c (Q). By the 
same proof of g2D, as A — >• A*, u\ — >■ u A . in n Lf+^fi). Then 

r^" n / (w A , + l) 2p < lim r^" n / (u A + l) 2p 



B r (as) A ^ A * ./B^x) 



< C lim r 4 ^"" / (u A + l) p+1 + (Au A ) 2 

JB 2r (x) 
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< Cr 4 ^- n f (u x * + l) p+1 + (Au x *) 2 . 

JB 2 r{x) 

Here we have used Fatou's lemma to deduce the first inequality. □ 
Below we denote u = u\* + 1. Inequality ( 16. 5ft implies that 

f u 2p <Cr n -^. (6.7) 

JB r {x) 

for any ball B r (x) C fl, with the constant C depending only on p and 
n. See for example the derivation of Eq.(2.16) in [32] • Similarly, u 
also satisfies (I4.3p for any ball B R (x) C Q. Estimate (I6.5P will play 
a crucial role in our proof of the e-regularity lemma. Note that both 
( 16. 5 P and ( 16. 7p are invariant under the scaling for ( II .ip . Theses two are 
also preserved under various limits (The precise notion of limit will be 
given below). 

To prove the partial regularity of u, first we need the following im- 
provement of decay estimate. 

Lemma 6.2. There exist two universal constants Eq > and 9 G (0, 1), 
such that if u is a positive stable solution of (II. ip satisfying the estimate 
( 16. 5p . and 

(2R)^~ n [ u p+1 + (Ah) 2 = e < e . 

Then 

(9R) 4 ^- n [ u p+1 + (Am) 2 < £ -. 

Proof. By rescaling, we can assume R—l. By ( 16.51) . we have 

f u 2p <C I u p+1 + (Ah) 2 < Ce. (6.8) 

■J B 3 / 2 J B 2 

By L 2 estimates applied to u, 

\\u\\ w ^(b 4/3 ) < c (lK|| L 2 (B3/2) + ||m||l2(b 3 /2)) < Ce ^- 

We can choose an r G (1,4/3) so that 

\\u\\w^2(dB ro ) < Ce^. (6.9) 
Now take the decomposition u = u\ + n 2 , where 

A 2 ui = u p , in B ro , 

u\ = Aui = 0, on dB ro (0), 
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and 



A 2 u 2 = 0, in B. 



u 2 = u, Au 2 = Au, on dB ro (0). 
By this decomposition, 



Am!Au 2 = 0. 



B rQ 



Hence 

(An) 2 = / (Am) 2 + / (Au 2 ) 



In particular, 

(Au 2 f < Ce. (6.10) 



Br 



By elliptic estimates for biharmonic functions and (16. 9p . we have 

1/2 



sup|u 2 |<C / u 2 + (Au) 2 ) <Ce— 

B1/2 \JdB ro 

Since Au 2 is harmonic, (Au 2 ) 2 is subharmonic in B ro . By the mean 
value inequality for subharmonic functions and (I6.10p . for any r G 

(0,r ), 

r 4 ^- n I {Au 2 f < r 4 ^r -" ! {Au 2 f < Cr^e. 

For ui, first by the Green function representation (cf. Section 4.2 in 
[18]), we have 



p+i 



IMU^) < c\\u p \\ LHBro) < C ^ u p+1 j < Ce-K (6.11; 

Then by L 2 estimates using (16. 7p . we have 

IKIk4, 2(Bi . o) < c {\\u p \\mB ro ) + IKIUi(s ro )) < Cei 

By Sobolev embedding theorem, we have 

||Mi|| 2n < C£2. 

Then an interpolation between L and L™- 8 gives 

IKIU^) <C^ +2<5 , 

where 5 > is a constant depending only the dimension n. 
Next, by interpolation between Sobolev spaces, we get 

||Aui|| La(Bn) ) < £~ 5 |k|U2 (Bro) + C , £ <5 ||A 2 n 1 || L2{Bro) < 
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Multiplying the equation of u\ by u\ and integrate by parts, we get 

.1+25 



/ u p Ul = I (Anx) 2 < Ce 1 

^ r ^ r 



By convexity, there exists a constant depending only on p such that 
u p+1 < C (u p+1 + u p 2 +1 ) . 
For r G (0, 1/2), which will be determined below, 

4£±l-n f p+1 ^ n 4£±i n. f p+1 . n 4£±i n f p+1 

r p- 1 / u p < Cr p- 1 / Ui + Cr p- 1 / 

J Bj- J By J Bj- 

< Cr 4 ^~ n I {u + u 2 y Ul + L7r 4 f^sup|M 2 r +1 

J B r B r 

< CrVi'" / u p Ui + CV p^i"™ / e^+Tm + Cr 

J By J By 

„ 4£±i-n /" p . ^ 4£±i-n /" +- , n 4^±i 

< Cr p- 1 / u p u\ + Cr p- 1 / £p +i Ui + Cr ?-*£ 

>/ _B ro J B ro 

„ 42±l-n 1+25 , /-y 4£±i-n . ^ 4^±i 

< Cr p- 1 £ + Cr p- 1 5P+ 1 +Cr v~ x £. 
For (Am) 2 , we have 

r^~ n ! (An) 2 < Cr 4 ^"" / (Am) 2 + Cr 4 ^"" / (Au 2 ) 2 

Cr 4 ^" n / (Aui) 2 + Cr 4 ^r "" / (An 2 ) 2 



< 



n 4£±i-n 1+25 , n 4£±i 

< Cr p- 1 e + Cr p- 1 £. 



Putting these two together, we get 

4£±i-n f /a \2 , p+1 / ^ 4£±i-n 1+25 , n 4£±i-n , /-y 4^±i 

r p- 1 / (Am) + u p < Cr p- 1 e T + Cr p- 1 5P+ 1 + Cr p-!£. 

J B r 

We first choose r = 9 e (0, 1/2) so that 

Cfl 4 ^ < 1. 
~~ 4 

Then choose an £ so that for every e e (0, £q), 

^4|±i-„ £ l+25 tf^-n^ < 1 £ 

~~ 4 

By this choice we finish the proof. □ 
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Remark 6.3. Lemma lKB also holds for a sign- changing solution u of 
(11.11) if it satisfies 

r^~™ f \u\ 2p < Cr^~ n [ \u\ p+1 + (An) 2 , (6.12) 

JB r (x) JB 2r (x) 

for any ball B2 r {x) C Q. For the proof, we need to introduce a new 
function u\, which satisfies 

A 2 mi = in B ro , 

U\ = Aui = 0, on OB ro (0), 

By the maximum principle, u\ > \u\\ > 0. By the same method for u\, 
we have 

( < Cs 1+2S . 

J B r() 

We can use this to control \u\ p \ui\. 

Lemma 6.4. There exists a universal constant e* > and 9 G (0, 1), 
such that if u is a stable solution of (11. ip satisfying (I6.12p . and 

(2i?) 4 ^" n I (Am) 2 + \u\ p+l = e<s , 

JB 2R (x ) 

then u is smooth in Br, and there exists a universal constant C(sq) 
such that 

sup |u| < C(s*)R~~. 

Br(x ) 

Proof. By choosing a smaller e*, we can apply Lemma [6.21 to any ball 
B r (x) with x G B R (x ) and r < -R/4, which says 

(6rf^- n [ (Am) 2 + \u\ p+1 < - r ^- n [ (Au) 2 + \u\ p+1 

J Bg r (x) 2 J B r (x) 

iterating the above implies 

/ (An) 2 + \u\ p+1 < Cr n ~^ +6 

JB r (x) 

for any x G B\ and r < 1/8. Here 5 > is a constant depending only on 
Eq and 9 in Lemma 16.21 In other words, u belongs to the homogeneous 

Morrey space L p+1 ' n+5 ~ i: p =I {Bi). Then the Morrey space estimate for 
biharmonic operator gives the claim, see the appendix. □ 

Together with a covering argument, this lemma gives a bound on the 
Hausdorff dimension of the singular set of u{= u x * + 1) 

dimS < n — 4^- — . 

p — 1 
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In particular, u is smooth on an open dense set. 

For any x G Q and A G (0, 1), define the blowing up sequence 

u x (x) = \p^u(x + Ax), A ->• 0, 

which is also a stable solution of (If .ip in the ball Bi/\(0). 
By rescaling (16. 7p . for all A G (0, 1) and balls B r (x) C -Bi/a, 

/ (u x ) 2p < Cr n ~^ . 

JB r (x) 

By elliptic estimates, u x is uniformly bounded in Wf^iW 1 ). Hence, up 
to a subsequence of A — > 0, we can assume that u x -> m° in W£ 2 (R n ) 
and Lf ^ 1 (R n ) (by the same proof of (14.41) ). By testing the equation 
for u x (or the stability condition for u x ) with smooth functions having 
compact support, and then taking the limit A — > +oo, we see u° is a 
stable solution of (O) in R n . 
We have 

Lemma 6.5. For any r > 0, E(r;0,u°) = lim E(r; 0, u). So u° is 

r->0 

homogeneous. 

Proof. A direct rescaling shows E(r;0,u x ) = E(Xr; Xq, u). By the 
monotonicity of E(r; Xq, u), we only need to show that, for every r > 0, 

E(r;0,u°) = lim E(r; 0, u x ). 

Because u x is uniformly bounded in W A ' 2 (B r ) and L 2p (B r ), by the com- 
pactness results in Sobolev embedding theorems and trace theorems, 
and interpolation between L q spaces (see (14.41) ). we have 

lim / (An A ) 2 = / (An ) 2 . 

A^+oo J Br J Br 

lim / (u x ) p+1 = f (u°) p+1 . 

A^+oo J Br J Br 

u x u° in W 2 ' 2 {dB r ). 

The last claim implies those boundary terms in E(r; 0, u x ) converges 
to the corresponding ones in E(r; 0, u°). Putting these together we get 
the convergence of E(r; 0, u x ). 

Since for any r > 0, E(r;0,u°) = const., by Corollary 12. 1[ u° is 
homogeneous. □ 

Here we note that since u satisfies (14.31) for any ball Br{x) C fi, so 
by the same argument as in the proof of Lemma I4.4[ we can prove that 
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E(r; x, u) is uniformly bounded for all x and r £ (0, 1). Since E(r; x, u) 
is non-decreasing in r, we can define the density function 

Q(x,u) := lim E(r; x, u). 

Lemma 6.6. (1) G(x,n) is upper semi- continuous in x; 

(2) for allx, Q{x,u) > 0; 

(3) x is a regular point of u if and only Q(x, u) = 0/ 

(4) there exists a universal constant e > 0, x £ S(u) if and only if 

G(x, u) > Sq. 

Proof. By the W^ 2 regularity of u, for any r > fixed, E(r; x, u) is 
continuous in x. Q(x,u) is the decreasing limit of these continuous 
functions, thus is upper semi-continuous in x. 

If u is smooth in a neighborhood of x, direct calculation shows 
Q(x, u) = 0. Since regular points form a dense set, the upper semi- 
continuity of G gives G > 0. 

By Lemma [6 A\ if x is a singular point, for any r > 0, 

/ (Au) 2 + u p+1 >e r n -^. 

JB r (x) 

In other words, for any A > 0, for the blowing up sequence u x at xo, 

[ (An A ) 2 + (u x Y +1 > e . 

JBi(p) 

Then because u x — > u° in W^(R n )nL^(R n ) (see the proof of Lemma 

ESD, 

/ (An ) 2 + (u°) p+1 = lim / (An A ) 2 + (n A ) p+1 

= limA- n+4 ^ f (An) 2 + {u) p+1 > e . 
Jb x (o) 

Hence u° is nontrivial, and by Remark 13.21 and Lemma 16. 5[ 

e(x,u) = E(l;0,u°) > c{n,p)e . 

Here c(n, p) is a constant depending only on p and n. 

On the other hand, if Q(x,u) < c(n,p)eo, then by Remark |3.2[ for 
any blow up limit u° at x, 

f (An ) 2 + (u°Y +1 < e . 
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Then by the convergence of u x in iy / 2 ' c 2 (M n )nLf+ 1 (R n ), for A sufficiently 
small, 

A 4 ^"" / (An) 2 + u p+1 = [ (Au x ) 2 + (u x ) p+1 < e . 
JB x (x) Jbi(o) 

By Lemma E31 u is smooth in B\/ 2 (x). Consequently, Q(x,u) = 0. 
These finish the proof of the last two claims. □ 

Remark 6.7. If lim u x = u° in some sense (for example, as in the 

A^O 

above blowing up sequence) so that for any x andr > 0, lim E(r; x, u x ) = 

A^O 

E(r; x, u°), then 

lim6(x,M A ) < e(x;u°). 

A— >0 

That is, Q(x; u) is also upper semi- continuous in u. 

Lemma 6.8. Let u £ Wf^iW 1 ) D L\^ c (R n ) be a homogeneous stable 
solution of (11.11) on W 1 , satisfying the monotonicity formula and the 
integral estimate (16.71) . then for any x ^ 0, Q(x,u) < 6(0, u). More- 
over, ifQ(x,u) = 0(0, u), u is translation invariant in the direction x, 
i.e. for all t £ R, 

u(tx + ■) = u(-) a.e. in R™. 

Proof. With the help of the integral estimate (16.71) . similar to Lemma 
S3, for any x £ R n , 

lim E(r; Xo, u) < C. 

And we can define the blowing down sequence with respect to the base 
point xo, 

u x (x) = \~u(x + Aa;) A — > +oo. 
Since u is homogeneous with respect to 0, 

u x (x) = n(A" 1 x + x), 

which converges to u(x) as A -> +oo in P^ 2,2 (R n ) n Lf+^R"). Then 
Lemma 16.51 can be applied to deduce that 



9(0; w) = E(l;0,u) = lim £(l;0,u A ) 

A— s-+oo 

= lim E(X;xo,u) 

A— >+oo 

> Q(x ;u). 

Moreover, if 6(xo; u) = 0(0, u), the above inequality become an equal- 
ity: 

lim E(X; xo, u) = O(xo; u). 

A->+oo 
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This then implies that E(X; Xq, u) = Q(xq] u) for all A > 0. By Corol- 
lary HHJ u is homogeneous with respect to xq. Then for all A > 0, 

u(x + x) = Ap T7I m(x + Ax) = «(A _1 x + x). 

By letting A — > +oo and noting that u(\~ 1 xo+-) are uniformly bounded 
in Wf G ' c 2 (M n ), we see 

u(x + •) = m(-) a.e. on lR n . 

Because u is homogeneous with respect to 0, a direct scaling shows 
that 6(txo; u) = 6(xo; u) for all t > 0, so the above equality still holds 
if we replace Xq by txo for any t > 0. A change of variable shows this 
also holds if t < 0. □ 

With this lemma in hand we can apply the Federer's dimension re- 
duction principle (cf. Appendix A in [30]) to deduce Theorem 11.91 

Appendix A. Proof of Estimate in Lemma 16.41 
Let us use the notation 

,, 7 ,o = sup ( r~ 7 / |/f 



x,r 



JB(x,r)nn 

L«"(fi) = \u G L 9 (fi) : ||u||, )7 ,n < oo}, 

where Q C M n is a bounded domain, 0<7<n, l<g< oo. 

For completeness we give a proof of the following result, which is an 
adaptation of [22j [23] . 

Lemma A.l. Assume u is a weak solution of 

A 2 u=\u\ p ~ 1 u inBi(0) 

and u e L p,n ^ 4 5^" H (5 1 (0)) /or some 5 > 0. T/ien it zs bounded in 
£1/2(0). 

We need some preliminaries. Let 



/«(/)(*)= / k - </P +c Y(y) 

Lemma A.2. (|221 Lemma 1]; If f E L 1,7 (R n ), < e < 7 and 1 < 
p < n ~ e then 

\I a {f)\ p {x)dx <Cdiam(tt) n - e - {n - a - e)p [ \f\dx (A.l) 
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Lemma A. 3. (Campanato pQ,) Let < 7 < n and c > 0. Assume 
(j) : (0, R] — > R is a nonnegative nondecreasing function such that 

<c(— 0(r) + r 7 ) for all < p < r < R. 



Then there is C depending only on n, 7, c such that 

0(p) < Cp 7 (^ + 1) for all < p < r < R. 

Lemma A. 4. Let v satisfy A 2 v = in 5^(0). Then there is C such 
that 

\v{x)\ < — I \v\dy for all \x\ < \r. (A.2) 
K Jb r (o) 1 

Proof. By scaling we can restrict to R = 1 and v G C 4 (5i(0)). Let 
77 G C7°°(]R n ) be a cut-off function with r](x) — 1 of \x\ > \ and 77(2;) = 
for \x\ < ~. Let T(x) = c n \x\ n ~ 4 be the fundamental solution of A 2 in 
W 1 , c n > 0. Then 

v(x) — I v(y)A 2 (T(x — y)r](x — y))dy for \x\ < - 

and (CO) follows. □ 

Proof of Lemma \A.1[ Let Ri < 1 (close to 1), |x| < i?i and < r < 
^-j^-. Let Mi = T * (lul^uxBrix)) where T(x) = c n |x| n ~ 4 is the funda- 
mental solution of A 2 in M. n , c n > 0, and XfirO) * s the indicator function 
of £ r (x). Let u 2 =u-u 1 . Then A 2 w 2 = in B r (x). By flA~2l) 

C f 

\u 2 (z)\ < — \u 2 \ for z G B r / 2 {x). 



B r (x) 



Let ?/ G B r / 2 (x) and < p < |. Integrating in B p (y) and using Holder's 
inequality 



Therefore 



\u 2 \P<C(^) n \u 2 \ p . 

B P (y) r JB r (x) 



I \u\ p <C [ |«i| p + C(~) n / \u 2 \ 

JB p {y) JB p (y) r J B r (x) 

1 

'B r (x) JB r {x) 



r 



< C(^) n / \u\ p + C / M p . (A.3) 
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Let 7q = n — 4^- + 5. Using OA.ip with a = 4, 7 = 7 and e a number 
such that n — 4-^- < e < 70 we have 



B r (ic) J B r (x) 



Then, combining with (I A . 3 1) we obtain 



u |p < C{-) n [ \u\ p + Cr n - £ - (n - 4 - e)p 



\ u \ 

B p (y) r ' JB r (x)' JB r (x) 

< (J(Py i / | U |P _)_ (J r n-£-(n-4-e)p+y 

r JB r (x) 

for any y G B r / 2 (x), < p < |. We have the validity of the inequality 
for < p < r, possibly increasing C. Using the Lemma of Campanato 
(Lemma IA.3j) . 

/ \u\ p < C p n - e -( n ~' l - € )p+'yo 

JB p {y) 

for < p < r, which shows that u G L pni (B^ ) where i?i < 1 can 
be chosen arbitrarily close to 1, and 71 = n — e — (n — 4 — e)p + 70 
can be chosen arbitrarily close to n(l — p) + 4p + 70. In particular we 
can choose 71 > 70- Repeating the process, we can find a decreasing 
sequence Ri — > g and an increasing sequence 7» — >• n — 4 such that 
u E L p ^(B Ri ). Then by Lemma \AJ\u E L q (B 3/A (0) for all q > 1 and 
by standard elliptic regularity m G L°°(Si/ 2 ). □ 



References 

[1] S. Campanato, Equazioni ellittiche del Ildeg ordine espazi 

Mat. Pura Appl. (4) 69 1965 321-381. 
[2] A. Chang, L. Wang, P. Yang, A regularity theory of biharmonic maps, 

Comm. Pure Appl. Math. (9) 52 1999 1113-1137 
[3] C. Cowan, P. Esposito and N. Ghoussoub, Regularity of extremal solu- 
tions in fourth order nonlinear eigevalue problems on general domains, 

DCDS-A 28 (2010), 1033-1050. 
[4] C. Cowan and N. Ghoussoub, Regularity of semi-stable solutions to 

fourth order nonlinear eigevalue problems on general domains, to appear 

Cal. Var. PDE DOI 10.1007/s00526-012-0582-4. 
[5] J. Davila, I. Flores and I. Guerra, Multiplicity of solutions for a fourth 

order equation with power-type nonlinearity, Math. Ann. 348 (2010), 

143-193. 

[6] Dancer, E. N. Moving plane methods for systems on half spaces. Math. 
Ann. 342 (2008), no. 2, 245-254. 



BIHARMONIC SUPERCRITICAL PROBLEMS 



39 



[7] J. Davila, L. Dupaigne, I. Guerra and M. Montenegro, Stable solutions 

for the bilaplacian with exponential nonlinearity, SI AM J. Math. Anal. 

39 (2007), 565-592. 
[8] J. Davila, L. Dupaigne and A. Farina, Partial regularity of finite Morse 

index solutions to the Lane-Emden equation, J. Fund. Anal. 261 (2011), 

218-232. 

[9] L. Dupaigne, Variations elliptiques. Habilitation a Diriger des 

Recherches. 11 december 2011. 
[10] L. Dupaigne, A. Farina and B. Sirakov, Regularity of the extremal 

solution, Geometric Partial Differential Equations , edited by Sc. Norm. 

Super. Pisa. To appear. 
[11] L. Dupaigne, M. Ghergu, O. Goubet and G. Warnault, The 

Gelfand problem for the biharmonic operator, Arch. Rat. Mech. 

DOI:10.1007/s00205-013-0613-028. To appear. 
[12] Fleming, Wendell H., On the oriented Plateau problem. Rend. Circ. 

Mat. Palermo (2) 11 1962 69-90. 
[13] A. Farina, On the classification of solutions of the Lane-Emden equation 

on unbounded domains of M. N , J. Math. Pures Appl. 87 (2007), 537-561. 
[14] F. Gazzola and H. -Ch. Grunau, Radial entire solutions for supercritical 

biharmonic equations, Math. Ann. 334 (2006), 905-936. 
[15] C. Gui, W. M. Ni and X. F. Wang, On the stability and instability of 

positive steady states of a semilinear heat equation in R . Comm. Pure 

Appl. Math. 45 (1992), 1153-1181. 
[16] Yuxia Guo, B. Li and J. Wei, Large energy entire solutions for the Yam- 

abe type problem of polyharmonic operator. J. Diff. Eqns. 254(2013), 

no.l, 199-228. 

[17] Z.M. Guo and J. Wei, Qualitative properties of entire radial solutions 
for a biharmonic equation with supcritical nonlinearity, Proc. American 
Math. Soc. 138 (2010), no.ll, 3957-3964. 

[18] F. Gazzola, H.-C. Grunau, G. Sweers, Polyharmonic boundary value 
problems. Positivity preserving and nonlinear higher order elliptic 
equations in bounded domains. Lecture Notes in Mathematics, 1991. 
Springer-Verlag, Berlin, 2010. 

[19] H. Hajlaoui, A. A. Harrabi and D. Ye, On stable solutions of biharmonic 
problem with polynomial growth. larXiv:1211.2223\ u2 (2012). 

[20] P. Karageorgis, Stability and intersection properties of solutions to the 
nonlinear biharmonic equation, Nonlinearity 22 (2009), 1653-1661. 

[21] Y. Li, C.-S. Lin, L. Nirenberg, Nonexistence results to cooperative sys- 
tems with supercritical expoenents in R™ , preprint. 

[22] F. Pacard, F. A note on the regularity of weak solutions of Au = u a in 
R n , n > 3. Houston J. Math. 18 (1992), no. 4, 621-632. 

[23] F. Pacard, Convergence and partial regularity for weak solutions of 
some nonlinear elliptic equation: the supercritical case. Ann. Inst. H. 
Poincar Anal. Non Linaire 11 (1994), no. 5, 537-551. 



40 



J. DAVILA, L. DUPAIGNE, K. WANG, AND J. WEI 



[24 
t 25 

[26 
[27 
[28 
[29 
[30 

[31 
[32 
[33 
[34 
[35 



F. Pacard, Partial regularity for weak solutions of a nonlinear elliptic 

equation. Manuscripta Math. 79 (1993), no. 2, 161-172. 

P. Polacik, P. Quittner, P. Souplet, Singularity and decay estimates in 

superlinear problems via Liouville-type theorems. I. Elliptic equations 

and systems. Duke Math. J. 139 (2007), no. 3, 555-579. 

P. Pucci and J. Serrin, A general variational identity, Indiana Univ. 

Math. J. 35(1986), 681-703. 

P. Pucci and J. Serrin, Critical exponents and critical dimensions for 
polyharmonic opera- tors, J. Math. Pures Appl. 69(1990), 55-83. 

F. Rellich, Perturbation theory of eigenvalue problems, Gordon and 
Breach Science Publisher, New York, (1969). 

G. V. Rozenblum, The distribution of the discrete spectrum for singular 
differential operators, Dokl. Akad. SSSR 202 (1972), 1012-1015. 

L. Simon, Lectures on geometric measure theory. Proceedings of the 
Centre for Mathematical Analysis, Australian National University, 3. 
Australian National University, Centre for Mathematical Analysis, Can- 
berra, 1983. 

B. Sirakov, Existence results and a priori bounds for higher order elliptic 
equations and systems, J. Math. Pures Appl. 89 (2008), 114-133. 
P. Souplet, The proof of the Lane-Emden conjecture in four space di- 
mensions, Adv. Math. 221 (2009), 1409-1427. 

J. Wei and X. Xu, Classification of solutions of high order conformally 
invariant equations, Math. Ann. 313(2) (1999), 207-228. 
K. Wang, Partial regularity of stable solutions to the supercritical equa- 
tions and its applications, Nonlinear Anal. 75 (2012), no. 13, 5238-5260. 
D. Ye, J. Wei, Liouville Theorems for finite Morse index solutions of 
Biharmonic problem, Math. Ann. to appear. 



J. Davila - Departamento de Ingenieri'a Matematica and CMM, Uni- 

VERSIDAD DE CHILE, CASILLA 170 CORREO 3, SANTIAGO, CHILE. 
E-mail address: jdavilaOdim.uchile . cl 

L. Dupaigne- LAMFA, UMR CNRS 7352, Universite Picardie Jules 
Verne, 33 rue St. leu, 80039 Amiens, France 
E-mail address: louis.dupaigne@matli.cnrs.fr 

K. Wang- Wuhan Institute of Physics and Mathematics, The Chinese 
Academy of Sciences, Wuhan 430071, China. 
E-mail address: wangkelei@wipm.ac.cn 

J. Wei - Department Of Mathematics, Chinese University Of Hong 
Kong, Shatin, Hong Kong and Department of Mathematics, University 
of British Columbia, Vancouver, B.C., Canada, V6T 1Z2. 

E-mail address: wei@math . cuhk . edu . hk 



